Electromagnetically induced transparency and reduced speeds for single photons in a 

fully quantized model 
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We introduce a simple model for electromagnetically induced transparency in which all fields 
are treated quantum mechanically. We study a system of three separated atoms at fixed positions 
in a one- dimensional multimode optical cavity. The first atom serves as the source for a single 
spontaneously emitted photon; the photon scatters from a three-level A-configuration atom which 
interacts with a single-mode field coupling two of the atomic levels; the third atom serves as a 
detector of the total transmitted field. We find an analytical solution for the quantum dynamics. 
From the quantum amplitude describing the excitation of the detector atom we extract information 
that provides exact single-photon analogs to wave delays predicted by semi-classical theories. We 
also find complementary information in the expectation value of the electric field intensity operator. 



I. INTRODUCTION 



Controlling the phase coherence in ensembles of mul- 
tilevel atoms has led to the observation of many strik- 
ing phenomena in the propagation of near-resonant light. 
These phenomena include coherent population trapping, 
lasing without inversion, electromagnetically induced 
transparency, and anomalously slow and anomalously 
fast pulse velocities. Resonant enhancement of the index 
of refraction without accompanying increased in absorp- 
tion was proposed and observed [||, ||] in 1991, and 
drastic reductions in the group velocity of pulses were 
discussed shortly thereafter Q . Recent experiments have 
taken the reduction of the speed of light to extreme limits 
H |(| and at the other extreme lie observations of seem- 
ingly superluminal light |?], ||. An overview of recent 
developments in the control of photons is presented in 
Ref . H . An earlier review of electromagnetically induced 
transparency was presented by Harris jl(| , while Lukin et 
al. present a overview of phase coherence in general, with 
an extensive list of references pL Such effects are also 
discussed in a recent texts (see, for example, Ref. Q). 
In most previous work the phenomenon of electromagnet- 
ically induced transparency and the accompanying dras- 
tic slowing of the speed of light are treated using semi- 
classical theory in which the atoms of the medium are 
treated quantum mechanically and the fields are treated 
classically. We use a model in which the entire system is 
treated quantum mechanically, and study the propaga- 
tion of a field state containing a single photon. Although 
coherent states of a single mode quantized field are of- 
ten considered as the "most classical," the single photon 
states that we study exhibit striking parallels with clas- 
sical fields. 



We study a system of three separated atoms at fixed 
positions in a one-dimensional multimode optical cavity. 
The first atom serves as the source for a single spon- 
taneously emitted photon; the photon scatters from a 
three-level A-configuration atom which interacts with an 
additional single-mode field coupling two of the atomic 
levels; and the third atom serves as a detector of the total 
transmitted field. In the spirit of Feynman's derivation 
of the classical index of refraction from the interaction of 
a field with a single oscillator Jl^ ], we infer the proper- 
ties of a medium exhibiting electromagnetically induced 
transparency from the interaction of the spontaneously 
emitted quantum field with the single quantized scatter- 
ing atom. We find an analytical solution for the quantum 
dynamics, including reradiation from the scatterer, and 
from this we deduce quantum delays that characterize 
the propagation of the field. These delays are equiva- 
lent to those predicted by semi-classical theories. In our 
quantum model all delays are clearly the result of inter- 
fering amplitudes that reshape the temporal envelope of 
the probability of detecting the transmitted photon. This 
effect is most clearly illustrated in the graphs of detection 
probability vs. time displayed in Sec. M. This work is an 
extension of the model we have used previously to study 
quantum manifestations of classical wave d elay s induced 
by scattering from simple two- level atoms jL4j . We note 
that the analytical results obtained in this paper may be 
verified using the straightforward numerical techniques 
like those used in Refs. |l|, |l8|. 
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II. REVIEW OF SEMI-CLASSICAL THEORY 

Electromagnetically induced transparency can be ob- 
served in the simple three-level atom illustrated in Fig. [I]. 
A strong laser "coupling" laser with angular frequency w c 
is tuned to resonance with the transition between atomic 
levels B and C, while a weak "probe" laser with angular 
frequency u> p excites the transition between levels A and 
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FIG. 1: Level scheme for simple electromagnetically induced 
transparency. A strong field resonantly couples levels B and 
C, and weak probe is near resonance with the transition be- 
tween levels A and C. The upper level C decays via sponta- 
neous emission at a rate 7 to the ground state A. All other 
damping mechanisms are assumed to be negligible. 



C. In this paper we consider the simplest case in which 
decay from level C to B is small enough to be neglected, 
and the only damping is due to emission at a rate 7 from 
level C to the ground state A. (In the this paper all decay 
rates 7^ refer to decay of probability of finding the atom 
in the excited state.) We also assume that there is no in- 
coherent pumping driving population between the levels 
of the atom. 

In the limit of a weak probe, the complex susceptibility 
is given by [pi 



X 



w£/4- S 2 -i6*f/2' 



(1) 



where N is the density of the atoms, d is the dipole mo- 
ment of the transition between levels A and C, lor is the 
Rabi frequency of the coupling transition between levels 
B and C, and the detuning of the probe laser frequency 
from resonance is 



S = lu p - U)AC- 



(2) 



The real and imaginary parts of the susceptibility are 
illustrated in Fig. @ for a coupling field strength such that 
= 7/2. The index of refraction and the absorption 
coefficient can be calculated from the real and imaginary 
parts of the complex susceptibility respectively. When 
the probe field is resonant with the transition between 
levels A and C, i.e., 5 = 0, the absorption goes to zero 
and the index of refraction is a rapidly varying function 
of probe frequency. 

For later comparison with the results of our quantum 
model, we consider a classical monochromatic plane wave 
of frequency to which is normally incident on a thin slab 
containing atoms with the level structure illustrated in 
Fig. |l} The plane of the slab is normal to the z axis, 
the thickness of the slab is Az, and the density of the 
atoms is N. If the incident field is E t = E e- %ul{t - z / c \ 
the transmitted field on the far side of the slab is 



E t = E t e 



iux^z/(2c) 



(3) 



and for weak scattering the transmitted field is approxi- 
mately 
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(4) 



The scattering is characterized by the dimensionless pa- 
rameter N Az\d\ 2 uj I (Tieocy) . 

The dispersion in the response of the atoms in this 
model leads to delays of pulses traversing such a slab. 
The delay in the arrival of the peak of a modulation en- 
velope of a quasi-monochromatic pulse that is away from 
a region of anomalous dispersion is determined by the 
group velocity v g = du/dk = c/(n + ujp-), and is given 
by 
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For detunings such that i5 <C ^>ac the group delay is 



At n = 



/ 2NAz\d\ 2 u AC \ 
\ he c J 

^ (u> R + 45 2 ) [(14 - 4S 2 ) 2 - 4 7 2 <5 2 ] 

[(w 2 ? -4 ( J2)2 + 4 7 2 (5 2]2 



(0) 



The functional form of the delay when lur = 7/2 is illus- 
trated in Fig. S. When the probe field is resonant with 
the the transition between levels A and C, not only is the 
absorption zero, but the group velocity can be extremely 
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FIG. 2: Real and imaginary parts of the complex suscepti- 
bility for the three-level system illustrated in Fig. [j]. The 
strength of the coupling field is such that lor = 7/4. 
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necessary to calculate explicitly the field £ in order to 
determine At £ 2 . For the specific fields considered in this 
paper we will show that the "temporal-center-of-gravity" 
delay happens to be equal to twice the group delay. 




FIG. 3: Group delay for a classical pulse in a medium with 
susceptibility given by Eq. (|l|). As in Fig. ^, the strength of 
the coupiing field is such that lor = 7/2. 



small, and this is manifested in the large positive group 
delay at S — in Fig. ||. The central peak in Fig. || be- 
comes taller and narrower as the strength of the coupling 
field is reduced. We note that for some values of the de- 
tuning S the group delay is negative, which corresponds 
to group velocities greater than the vacuum speed of light 
c. Such "superluminal" velocities do not violate causal- 
ity, and are an effect of pulse reshaping by the dispersive 
medium. Similar reshaping effects and group velocities 
greater than c occur near simple two-level resonances in 
both classical and quantum theories |l4| . 

For pulses that are either not sufficiently monochro- 
matic, or not sufficiently far away from a region of anoma- 
lous dispersion, the simple concepts of phase and group 
velocity are inadequate to characterize all of the effects 
of pulse-reshaping as a field propagates. Several other 
velocities and delays have been developed (see, for ex- 



ample, Refs. ]Tj 



20 ) and in this paper we focus on a 
temporal center of gravity" of 



delay determined by the 
the field intensity of a pulse at a fixed position z "down- 
stream" from the slab containing the atoms comprising 
the medium, i.e., 
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/ t£(z,t) 2 dt\ 
J £{z 1 t) 2 dt J a ft er medium 
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(7) 



This is closely related to concepts used to define the cen- 
trovelocity in We have investigated this delay in 

classical and quantum mechanical models of scattering 
from simple two-level atoms in a previous paper |l4j . For 
quasi-monochromatic pulses far from resonance this de- 
lay is equivalent to the group delay, but in general it is 



III. QUANTUM MECHANICAL MODEL 

The quantum mechanical system we consider is illus- 
trated in Fig. [|, and consists of three atoms in a multi- 
mode one-dimensional optical cavity that extends from 
z = to z = L. This multimode cavity is oriented hor- 
izontally in the schematic representation of Fig. || The 
middle atom has an additional interaction with a single- 
mode field contained in the vertical cavity. The field in 
the multimode (horizontal) cavity plays the role of the 
probe field, and the field in the single mode (vertical) 
cavity represents the coupling field. (The finite optical 
cavities do not contribute to the physical phenomena un- 
der investigation; they simply provide a convenient quan- 
tization volume for the field modes used in our calcula- 
tion.) In the remainder of this section we will discuss 
the details of our model and the standard quantum op- 
tical Hamiltonian we use. We also present the analytical 
solution for the time dependence of the system. 

The atom on the left (atom 1) is a two-level atom which 
is initially in the excited state, and will be the source of 
the probe field. The middle atom (atom 2) which will 
scatter the radiation emitted by the source is a three-level 
atom with the "lambda" configuration of Fig. [j], with the 
energy difference between levels A and C close to that of 
the level separation of atom 1. (The highest level of all 
three atoms will be labeled as C.) Levels B and C of atom 
2 will interact with the single-mode coupling field which 
is assumed to be exactly on resonance. The coupling field 
will initially be in a state with a well-defined number of 
photons such that the Rabi frequency of the transition 
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FIG. 4: Quantum mechanical model consisting of three atoms 
at fixed positions. All three atoms interact with the radiation 
in the probe field in the multimode cavity (represented hori- 
zontally) and the middle atom also interacts with the single 
mode coupling field (represented vertically). 
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1 2 3 

FIG. 5: Level scheme of the source, scattering, and detector 
atoms. The source and detector are two-level atoms with 
identical transition frequencies. The scattering atom is a 
three-level atom in the A configuration. The transition fre- 
quency of the source and detector are detuned from the A-C 
transition of the scattering atom by an amount 5. 

between levels B and C is appropriate for the observation 
of electromagnetically induced transparency. The two- 
level atom on the right (atom 3) will serve as a detector. 
The detector atom is assumed to have the same resonant 
frequency as the source atom. The relative energy levels 
of all three atoms are illustrated in Fig. ||. 

We assume that the coupling strengths of the atoms 
to the probe field are all independent, so that the atoms 
have separate spontaneous decay rates. Although we find 
a solution for the dynamics for any decay rates, we focus 
our attention in this paper on cases in which the decay 
rate of atom 1 is very much smaller than the decay rate of 
atom 2. This condition insures that the spectrum of the 
radiation emitted by atom 1 will be very much narrower 
than the linewidth of the the scattering atom, and this 
allows us to compare our results from the quantum case 
with those of the semi-classical model which assumes that 
the sample is driven by a monochromatic wave. We also 
assume that the decay rate of the detector atom (atom 
3) is very much greater than any other dynamical rates 
in the problem. In this limit the excitation of atom 3 will 
closely follow the field drives it. For simplicity we also 
assume that the spontaneous decay rate from level C to 
level B (for atom 2) is negligible. 



The zero-field resonance frequencies of the A-C transi- 
tions of the atoms are labeled , where j — 1, 2, or 3, 
and the positions of the atoms will be labeled Zj . In the 
remainder of the paper we will assume that the atoms 
are at positions z\ — L/4, z 2 — L/2, and z 3 = 3L/4, as 
illustrated in Fig. ||, although our results for delay times 
do not depend on the exact positions. The standing wave 
field modes of the probe field cavity are separated in an- 
gular frequency by the fundamental frequency 

A = tt|. (8) 

This mode spacing may be small enough that many 
modes fall within the natural line-width of the atoms. 

For convenience we assume that the frequency of one 
of the modes corresponds exactly to the resonance fre- 
quency of atom 1 , the emitting atom, and that the length 
of the cavity is such that it contains an even number of 
wavelengths of this mode. We label the frequency of this 
mode ujq — fcoA, where fco is an integer divisible by 4. 
(This assumption affects the details of some of our cal- 
culations, but not our results concerning delay times.) 
The other mode frequencies will be enumerated from this 
mode so that 

u k = (k + k)A, (9) 

where k = 0,±1,±2, . . . 

As in the semi-classical case we wish to study the ef- 
fects of the detuning of the source field on the scattering 
of the radiation. We use the same symbol 5 as in the 
semi-classical case to represent the detuning of the field, 
but in the quantum case the detuning is directly tied to 
the properties of the source and scattering atoms: 

5 = ^ at) -4 at) - (10) 

Because the detector atom is assumed to have the same 
resonance frequency as the source atom we have = 

(at) 
^3 ■ 

We use as basis states the eigenstates of the atomic 
plus free-field Hamiltonian 



Hq — Hatoms + field 

— -f^atoms ~t~ -^probc ~t~ -^coupling 

= 7k^ at) |Ci>(<?i| + huji 2 at) \C 2 )(C 2 \ + h{Jf ] -u c )\B 2 )(B 2 \ +hLu { 3 at) \C 3 )(C a \+J2^kala k + hu c ala Cl (11) 

k 

I 



where a k and a\ are the lowering and raising operators 
for the fc th mode of the probe field, and a c and a\ act 
similarly on the single mode of the coupling field. (We 



have re-zeroed the energy scale to remove zero-point en- 
ergy of the field modes.) The basis states will be denoted 
as follows: 



■5 



• \C, A, A; 0, TV) — Atom 1 excited, atoms 2 and 3 in 
the ground state, no photons in the cavity modes, 
N photons in the coupling mode; 

• | A, C, A; 0, N) — Atom 2 in state C, atoms 1 and 3 
in the ground state, no photons in the probe field, 
N photons in the coupling mode; 

• \A, A, C; 0, N) — Atom 3 excited, atoms 1 and 2 in 
the ground state, no photons in the probe field, N 
photons in the coupling mode; 

• \A,A,A;lk,N) — All atoms in the ground state, 
one photon in the probe field mode with frequency 
(ko + k)A, N photons in the coupling mode; 

• \A, B, A; 0, N + 1) — Atom 2 in state B, atoms 1 
and 3 in the ground state, N + 1 photons in the 
coupling field. 

Each atom is coupled to all of the standing wave modes 
of the probe field, and we use the symbol gjk to label the 
coefficients characterizing the strength of the coupling of 
the j th atom to the k th mode of the probe field. In ad- 
dition atom 2 interacts with the coupling field, and this 
interaction is characterized by the constant g c . We use 
the standard electric-dipole and rotating-wave approxi- 
mation approximations to give the following interaction 
Hamiltonian: 

3 

^ int = S S h {9jka\\A 3 ){Cj | + g* jk a k \Cj) (A,-|) 

3=1 k 

+h (g c at\B 2 )(C 2 \ + n* c a c \C 2 )(B 2 \) . (12) 

The Rabi frequency lor of the transition between levels 
B and C is determined by the number of photons in the 
coupling mode and the coupling constant g c as follows: 



lor = 2g c VN + 1. 



(13) 



For convenience we assume that lor is real. 

We assume that the frequencies of all atomic transi- 
tions are very much greater than the fundamental fre- 
quency of the cavity, that is lo^ 3> A for all atoms, and 
similarly for lo c . In this limit we can make the approxi- 
mation that all modes that influence the dynamics of the 
system are near the atomic resonances, and the atom- 
field coupling constants can be factored into a product of 
a frequency-independent constant and a space-dependent 
coupling factor. The coupling constants gjk are given in 
terms of the electric dipole matrix element dj between 
the levels A and C of atom j, the effective volume of the 
cavity, V, and the permittivity of free space, eo, by 
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gjk 



d 



2he V 
Ctj sin [(ko + k)irzj/L] , 



sin [(ko + k)iTZj/L] 
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FIG. 6: Magnitude of the amplitudes for the atoms to 
be in the excited state, starting from the state |^(0)) = 
|e,<?,g;0, N). The decay rates of the atoms are 71 = 4,72 = 
64, and 73 = 1024; the Rabi frequency of the B-C transition 
of atom 2 is lor = 72/2 = 32, and 5 = 0. 



where in the last line we have defined the quantity 



Clj — dj 



LO 



(at) 



2he V 



(15) 



(14) 



which is independent of the cavity mode-frequency 

We write the state of the system as the linear combi- 
nation 

|V(t)> = c l (t)\C,A,A-$,N) + c 2 (t)\A,C,A-$,N) 
+c 3 (t)\A,A,C;<D,N) 
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+d(t)\A,B,A;<D,N+l) 

+ Y J b k {t)\A 1 A,A-l k ,N). (16) 

k 

Choosing the zero of the energy scale at the level of 
(uncoupled) state \C, A, A; 0, N), the Schrodinger yields 
the following set of coupled differential equations: 

ci = -~i^g lk b k , (17) 

k 

C2 = -i 92kbk + u B d/2 - Sc-^j , (18) 

c-3 = -i^293kb k , (19) 

fc 

d = -i(u) R c 2 /2-6d), (20) 

h = -i (y^g^Cj+kAb^j . (21) 

We solve this set of equations with the Laplace trans- 
form technique used by Stey and Gibberd Q for a 
Hamiltonian similar to that for a single atom at the 
center of one-dimensional cavity. Laplace transforms 
have also been used to solve the Schrodinger equation 
in similar problems with two interacting atoms in three- 
dimensions p2| , and in our own recent work on scat- 
tering from two-level atoms [[l4|. Because the Laplace 
transform technique is not new, and because we would 
like to focus on analogies with the fields in the semi- 
classical model and physical interpretation, we leave the 
details of our solution to the appendix, and simply quote 
our results here. 

The general features of the solution giving the time- 
dependences of the atomic excitation amplitudes are il- 



lustrated in Fig. |6j. The initially excited atom decays 
exponentially until the time t — 0.5L/c at which scat- 
tered and reflected radiation first returns to the atom. 
The amplitudes to find the other atoms excited are iden- 
tically zero until radiation first reaches them: the scat- 
tering atom first becomes excited at t — 0.25L/c and the 
detector atom is first excited at t = 0.5L/c. The three de- 
cay constants which characterize the spontaneous emis- 
sion rates of level C in each of the atoms emerge naturally 
in terms of the parameters of the Hamiltonian as 

7* = ^ = !^. (22) 
A c 

The causal nature of the dynamics is evident in that all 
disturbances are propagated at the speed of light c via 
the quantum field. The abrupt changes in the complex 
amplitudes at intervals of 0.5L/c are a manifestation of 
the finite speed of light and the atomic spacing of 0.25L. 

The abrupt changes appear in our analytic solution for 
the complex amplitudes Cj(t), d(t), and bk{t) as sums of 
terms with step functions that "turn on" at successively 
later intervals of 0.5L/c. In the following formulas giv- 
ing these amplitudes we truncate the expressions so that 
only the first excitation of atoms 2 and 3 are included. 
We also note that the following equations are specific in 
some details to the atomic positions Z{ used in this pa- 
per. The "turn-on" times and relative phases of terms 
will change with different positions, but the conclusions 
of this paper concerning delay times are unaffected by 
these details. The time dependence of the system is given 
by the following set of amplitudes, in which we use the 
labels Ti as a shorthand for complicated factors that are 
functions of ji, S, lur, and (in the case of bk), k: 



Ty + ( t - - ) F 2 



+e 



+ ■ 



ci(i) 
<*(*) 

C3(*) 



d(t) = e(t- ^ { e -¥(*-i)^ u + e KWa^ 

h{t) = (e-^M - e-*^) g lk T lb + Q (t - g 2k [e^^F^ + e-MM)^ 



[-(72-VTf-4<4)/4+i5](t-!) jP 3 + e [-(72 + V7 2 2 -4o ) | ? )/4+ i A-](t-i) ^ j 



e 1 1- 

+e -T( t -2)jT 11 L + ... 



(23) 
(24) 

(25) 
•(26) 
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+e [-(72-VT|-4a.|)/4+i5](t-i) ^ + e [-(72 + V7 2 2 -4o;|)/4+i«](t-^) p\ + 



(27) 



r 



Complete expressions for the factors Ti may be found in 
the appendix. 

Four time scales are manifest in the exponential fac- 
tors in these these equations: the exponential decay rate 
of the source atom 71 , which is assumed to be small, the 
exponential decay rate of the detector atom 73, which 
is assumed to be large, and two decay rates associ- 
ated with the scattering atom, ^72 + a/tJ — ^r) /2 an d 

(72 - \/72 - 4w |) A whcn thc 

coupling field strength 

is small, the last of these rates approaches zero. This 
means that decay from the scattering atom can be very 
slow, and this slow reradiation is at the origin of the 
reduced velocities which characterize media with electro- 
magnetically induced transparency. 

In the following sections we will focus on two quanti- 
ties: cs(t), the amplitude to find the detector atom ex- 
cited, and (£ 2 ) the expectation value of the square of the 
electric field operator, which is proportional to the field 
intensity. (The expectation value of the field operator 
itself is for any single-photon state.) In our investi- 



J 



gation of c 3 (t) we will consider only the first term in 
Eq. pa) describing the initial excitation of the detector 



atom. Similarly, we will investigate (£-) in regions to 
the right the scattering atom, and at times that exclude 
multiple scattering effects. 

It is useful to rewrite c 3 (t) as the sum of two pieces, the 
amplitude c° (t) for atom 3 to be excited in the absence of 
the scattering atom (or, equivalently, when 72 = 0), plus 
c|(t), the amplitude that is attributable to scattering: 



c 3 (t) = c 3 (t)+cl(t). 



(28) 



Setting 72 = in Eq. (Oa) gives 



4® 



e t 



1 
2 



V7i73 
7i - 73 



-*(t-i) 



(29) 



and subtracting this from Eq. (p5|) gives 



4(t) = e(*-i) { e -^(M).F2 + e [-(^^F^)/^^ 

+e -¥(M)^- 23 j (30) 



where 



^2. 



^21 = 



•7"22 



V7173 

7i - 73 J 7i - 7i 72 - 4<5 2 + uj 2 r + i2S(2ji - 72) 



72(71 + i2S) 



72 (72 - VtI - 4w ff) 



2^7173 \ 

. VH ~ (271 ~ 72 + y/H ~ 4c4 + *4j) (273 - 72 + ^7 2 2 ~ + 

2^777? \ 72 (72 + yg - 4^) 

. VH ~ H/ (271 - 72 - V^F^r + iiS ) (273 - 72 - V72 2 - + *4<5) 

V7i73 \ 72(73 +i2S) 

71-73 / 73 - 7273 - 4<5 2 + wfj + i2S(2-f 3 - 72) 



(31) 
(32) 

(33) 
(34) 



In our investigations of the quantum field itself we use the electric field operator in the form given by Meystre and 
Sargent p4| , and write the expectation value of the square of the field as 



(£ 2 ) = m)\ \ E («* + 4) [(/co + k)^] \ \m) 



(35) 



In the limit considered in this paper we can replace the frequencies uj^ under the radical with the constant 
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(at) 



After expanding the state vector as in Eq. (fL6|) and 



evaluating the sums, the expectation value can be written 
in terms of the amplitudes b k (t) to find the photon in the 
various cavity modes: 



(£ 2 




^ b k (t) sin (k + k) 



L 



(36) 



(In this expression we have dropped the infinite term 
arising from vacuum expectation value off 2 .) Evaluation 
of (£ 2 ) gives a space- and time-dependent representation 
of the localization of the energy of the photon Jl4], [l^, ^g) . 

The expression for (£ 2 ) in Eq. ( |36] ) is the square of a 
complex number that is analogous to the complex ana- 
lytic signal describing the classical field. We label this 
quantity £ q . m , i.e., 



f — 



(at) 



e V 



^b k (t) sin (fc, 



(37) 



We note that the quantity £ q . m . is effectively equivalent 
to (0\E( + )\ip 7 ) , where is the positive frequency part 



of the field quantized in terms of traveling wave modes, 
which has been identified as "the 'electric field' associ- 
ated with [a] single photon state" by Scully and Zubairy 
|l2| . In a previous paper fl4f| we extended the quantum- 
classical correspondence embodied in this quantity to 
fields that include scattering from two-level atoms, and 
here we extend it to include scattering from three-level 
atoms that make up a medium that exhibits electromag- 
netically induced transparency. 

For ease of comparison with previous results for the 
detector atom, we calculate the field at the fixed position 
3L/4 in a cavity that does not contain the detector 



atom. With no scattering atom present we find 1 14 




* (at) 

2e V 



(38) 



The energy density passing the point z = 3L/A exhibits 
an abrupt turn-on (because of the initial conditions we 
have chosen) followed by exponential decay [^2[ [j~4| [jl], 
p"6| . With a three- level scattering atom present at Z = 
L/2 we find 



-ie t 



72 72 



+ 



(- 




-(*-§) 



72(71 + i25)e i (* 2) 



7i 2 - 7i72 - 45 2 + uj 2 r + «25(27i - 72 ) 



Auj r 1 ' 



[-h2-v / l!-4u,%)/4+iS](t-±) 



V72 - ^R ( 2 7i - 72 + y/H - ^ R + M&) 
72 (72 + y/H - ) e h(7 2 +V^F^)/^](*-|) 



(271 



72 - vtF" 



4-1 



(39) 



r 



The first term in curly brackets is just the previous re- 
sult with no scattering atom present; the effect of the 
scattering is contained in the remaining terms. The 
steady state scattering is contained in the second term 
in curly brackets. This term decays at the slow rate 
71 reflecting the envelope of the incident radiation. We 
claimed previously that for large 73 the detector atom re- 
sponse reflects the field incident upon it, and comparison 
of Eq. (|3^), which gives the probability amplitude due 
to scattering for the detector atom, with Eq. (^9|), justi- 
fies this claim. In the limit 73 2> 72 , 71 the factors ^20 
through ^"23 that arise in Eq. ( p0| ) are proportional to the 
the coefficients in front of the corresponding exponential 
terms in Eq. (|39|). 

The quantum pulse of Eq. (p8[) has the classical analog 



4° L = eft-^ ce-^-K-iW-i), 



(40) 



where C is a constant. If such classical pulses are incident 
on a thin slab of material with thickness Az, density 
TV, and classical susceptibility x given by Eq. ([!]), then 
the classical transmitted field £ c \. is identical in form to 
Eq. (|39|), except that the terms due to scattering (i.e., 
those proportional to 72) have a magnitude characterized 
by the small dimensionless factor N Az\d\ 2 u! / \he0cj2) ■ A 
derivation of this result is included in an appendix. 



IV. COMPARISON OF SEMI-CLASSICAL AND 
QUANTUM MECHANICAL SCATTERING 

In the semi-classical model the result of weak scat- 
tering of a monochromatic field is contained in Eq. (Q). 
In the limit of narrow bandwidth probe (71 <ti 72) and 
rapid detector atom response (73 3> 71,72) our quan- 
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turn probability for the detector atom to be excited by 
the scattered field is approximated by the first term of 
Eq. ©, 



c s 3 (t) 



eu-;| e -T(*-i)jF 20 



1812/2 



4(t) 



u 2 R /4 - 5 2 - i8 l2 /2 



(41) 



This is exactly the same functional form as the term 
due to scattering in the formula for the classical field, 
Eq. ([I]). In the classical formula the dimensionless fac- 
tor NAz\d\ 2 uj/(heocj) is assumed to be small. In our 
quantum model the magnitude of the scattering is de- 
termined by Q 2 ( or equivalently 72 ) , which characterizes 
the coupling of atom 2 to the probe field. In our one- 
dimensional model the coupling to the incident field and 
the decay rate of atom 2 are both completely determined 
by the single parameter , which means that it is not 
possible to make the effect of the scattering small with- 
out simultaneously making the linewidth of atom 2 very 
narrow. In a fully three-dimensional model the decay 
rate of atom 2 would be the result of the atom's coupling 
to many more modes, and not just those containing the 
incident field. This scattering into other modes would re- 
duce the scattering in the forward direction (the direction 
of the detector) from the amount predicted in our simple 
model. If the forward scattering is reduced by a factor /, 
then a more realistic expression for the excitation of the 
detector atom is 



c 3 (t) = c° 3 (t) + fc s 3 (t)- 



(42) 



In this formula / plays a role analogous to the dimen- 
sionless factor N Az\d\ 2 cu / (heocy) in the formulas for the 
transmitted field derived from semi-classical theory. 



V. TEMPORAL-CENTER-OF-GRAVITY DELAY 

The group velocity of a classical field has a clear inter- 
pretation for quasi-monochromatic pulses whose central 
frequency is far from a region of anomalous dispersion: 
it is the speed at which the peak of the modulation enve- 
lope travels. Group delays refer to the delay in the arrival 
of the peak of a pulse compared to the time expected 
for propagation through a vacuum. The pulses investi- 
gated in this paper have sharp leading edges, and this 
lack of a smooth modulation envelope means that the re- 
sults of simple classical theory for quasi-monochromatic 
pulses should not be expected to be a sufficient guide 
to full understanding. In this section we will investi- 
gate "temporal-center-of-gravity" delays in several clas- 
sical and quantum mechanical quantities. 



The first delay we investigate is derived from c 3 (t), 
the amplitude for the detector atom to be excited. As 
we have argued previously, this amplitude reflects the 
strength of the incident field in the limit that the response 
time of this atom is very short compared with other time 
scales, i.e., 73 3> 71,72- The effect of the scattering on 
this amplitude is evident in Fig. ^ in which | C3 (t) | 2 , the 
probability for the detector atom to be excited, is plotted 
for three values of the detuning 8, and also for the case 
in which no scattering atom is present. (The results are 
plotted for the specific case in which lur = 72/2.) For all 
detunings, C3(t) = for all times earlier than t = 0.5L/c, 
as is expected; all effects on the detector atom occur at 
times that preserve causality. The qualitative shapes of 
the detector response depend critically on the detuning 
8. For 8 = the steady-state absorption coefficient is 
zero, and this is reflected in the fact that at large times 
there is little effect of the scattering on the probability. 
At shorter times the effect of transient oscillations have 
a pronounced affect, and he temporal center of gravity of 
the response is clearly shifted toward later times relative 
to that of the response with no scattering atom present. 
The temporal center of gravity shifts to earlier times as 
the detuning is increased, and for values in the vicinity 
8 = 72/2 it occurs earlier than the in the case of no scat- 
tering. This is the region in which negative values of the 
group delay occur, as is illustrated in Fig. ||. The finite 
response time of the medium results in large brief trans- 
mission of leading edge of the pulse of the field before the 
high attenuation of the steady state sets in. (A similar 
effect occurs in scattering from simple two-level atoms 
fl4j| .) At larger values of the detuning, such as the case 
8 = 72 illustrated in the bottom graph of Fig. [?], most of 
the probability is again removed at early times, shifting 
the temporal center of mass back to later times. 

We quantify the ideas illustrated in Fig. by iden- 
tifying an effective arrival time of the photon with the 
temporal center of gravity of the probability that the de- 
tector atom is excited, i.e., 



rival 



Jt\c 3 (t)\ 2 dt 
J\c 3 (t)\ 2 df 



(43) 



(In evaluating the integrals in this equation, we use only 
the first term in the series of Eq. (E3) , and assume that 
the decay rates and distances are such that the effect of 
multiple scattering is negligible.) The delay imposed by 
the medium is then just the difference in the arrival times 
calculated with and without a scattering atom present, 



Ai C3 = 



ft\c 3 (t)\ 2 dt ft\4(t)\ 2 dt 
f\c 3 (t)\ 2 dt f\4(t)\ 2 dt- 



(44) 



To explore the effect of weak scattering we rewrite c 3 (t) 
in the form of Eq. (f42f), and assume that / < 1. Our 
quantum mechanical delay becomes, to first order in /, 



At,., = 



Jt\4(t) + fcm 2 dt Jt\c°(t)\ 2 dt 
J\c°(t) + fcm\ 2 dt J\4(t)\ 2 dt 
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FIG. 7: Probability for the detector atom to be excited as 
a function of time for various values of the detuning S. (All 
effects due to multiple scattering and reflections have been 
suppressed.) The effect of the "pulse reshaping" is dependent 
on the detuning. The temporal center of gravity of the prob- 
ability is shifted relative to that in the case of no scattering 
from relatively late times in the top graph to earlier times in 
the middle two graphs, and then back to later times in the bot- 
tom graph. The decay rates of the atoms are 71 = 4,72 = 64, 
and 73 = 1024; the Rabi frequency of the B-C transition of 
atom 2 is lur — 72/2 = 32. 



2/ 



JtRe[c° 3 (t)cm*] dt 
I\4(t)\ 2 dt 

ft\4(t)\*dtfRe[4(t)cm*} dt 

(I\4(tWdt) 2 



(45) 



It is straightforward (but involved) to evaluate the in- 
tegrals in Eq. ( |45|) using the expressions for c®(t) and 
c|(t) from Eqs. (|29|) and (|30"|). After taking the limit 
73 — > 00 and then letting 71 — > we arrive at the follow- 
ing expression for the quantum time delay: 



At r 



4/7: 



(lo 2 r + 4S 2 ) [(c4 - AS 2 ) 2 - 4 7 2 <5 2 



4S 2 ) 2 +4j 2 5 2 Y 



(46) 



Recalling that the parameter / characterizes the mag- 
nitude of the scattering in the forward direction in the 
same way that the quantity N Az\d\ 2 cu / (hegc^f) does in 
the semi-classical case, we see that our "temporal-center- 
of-gravity" delay time for this specific pulse is identical in 
functional form to the expression for the group delay of 
the classical field given in Eq. (g) . The magnitude of the 
temporal-center-of-gravity, however, is twice that given 
by the group delay, i.e., 



At C3 = 2(At g ) 



g ^classical • 



(47) 



The difference in the classical group delay and the 
temporal-center-of-mass delay should not be surprising, 
because the pulses studied in this paper do not satisfy the 
quasimonochromatic condition. The delay in the arrival 
time that we have defined is the result of the reshaping 
of the "pulse" of excitation of the detector atom. The ef- 
fect of scattering in the region of resonance is effectively 
to redistribute (in time) the probability that the detec- 
tor atom will be excited. The atom can rapidly absorb 
energy from the probe field, and the coupling field cre- 
ates a state of the system which can then return energy 
to the probe field at a rate which can be adjusted, via 
the Rabi frequency wr, to be arbitrarily slow. Although 
we have not demonstrated it explicitly in this work, we 
are confident that delays of quasi-monochromatic quan- 
tum pulses with smooth envelopes can be explained in 
the same manner. 

Because £ c \., £ q . m ., and cz(t) all have the same func- 
tional form (in the large 73 limit), it is easy to see that 
equivalent delays can be derived from the classical field 
using Eq. (R]), or from the quantum field using the anal- 
ogous equation 



At 



(s 2 



( ft{S(z = 3L/4y 



dt 



\ f{£(z = 3L/A) 2 )dt 

ht{£{z = ?,L/A) 2 )dt 
\j{£(z = 3L/4y)dt f 



with scatterer 



.(48) 



no scatterer 



VI. CONCLUSION 

We have considered a fully quantized model of scatter- 
ing by three-level atoms which can exhibit elect romagnet- 
ically induced transparency. Our model is simple enough 
that we are able to find an analytical solution describ- 
ing the complete dynamics of the system. Using this 
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model we have investigated propagation of single spon- 
taneously emitted photons through a medium exhibiting 
electromagnetically induced transparency, and have de- 
fined an effective time of arrival at a detector atom. We 
have also calculated the expectation value of the field 
intensity operator, and identified a quantum analog to 
the complex analytic signal describing the classical field. 
These quantum mechanical quantities exhibit delays that 
are clearly the result of pulse reshaping effects. We have 
compared our delays for quantized one-photon fields with 
those calculated for classical fields using the index of re- 
fraction determined from semi-classical theory, and find 
them to be in agreement. 



APPENDIX A: SOLUTION USING LAPLACE 
TRANSFORMS 

Before solving the coupled differential equations re- 
sulting from the Schrodinger equation, Eqs. ( |l7| - |2l] ) , we 
re-zero the energy scale and measure all energies rela- 
tive to the energy of the initial state |e,<?,<?;0) (ignoring 
^interaction)' Taking the Laplace transform of this set 
of coupled differential equations gives a set of coupled 
algebraic equations, 

i(sci(s)-l) = ^gikbkis), (Al) 

k 

{is + 6)£ 2 (s) = ^2g 2 kb k (s)+ujRd{s)/2, (A2) 

k 

isc 3 (s) = ^2g3kh(s), (A3) 

k 

isd(s) = w fi c 2 (s)/2- 5d(s), (A4) 
isb k (s) = ^ fl * fc c,(s)+fcA6 fc (s). (A5) 

3 

We solve these algebraic equations for the quantities 
Cj(s), d(s), and &fc(s), and then perform an inverse 
Laplace transform to recover the time dependence of 



Cj(t), d(t), and bk(t). The details of carrying out such 
calculations are quite involved, and were completed with 
the aid of a computer algebra system p5|. In this ap- 
pendix we outline our approach and present some of our 
intermediate results. 



We begin by solving Eq. (A4) for d(s) and Eq. ( A5) for 
6fc(s), and substitute the results in the first three equa- 
tions, giving 



sci - 1 = -iA (/nci + /i 2 c 2 + A1C3) 



f'2 



iA (A2C1 + f 22 c 2 + /23C3) 
4(s- iS) 



4(s- iS) 2 +uj 2 



R 



(A6) 



(A7) 



sc 3 = -iA (/13S1 + /23C2 + /33C3) , (A8) 
in which wc have defined the dimcnsionlcss sums 



flm - A2 



\_ \ ^ 9lkg* mk 
k £A 



(A9) 



In the limit in which the atomic resonance frequencies 
are very much greater than the fundamental frequency 
of the cavity, i.e., 3> A, these sums may be ap- 

proximated by extending the range for k from — 00 to 
+00, in which case the sums have relatively simple rep- 
resentations in terms of trigonometric functions. Explicit 
expressions for these sums are given in the Appendix of 
Ref. lO. 



After solving Eqs. (A6-A8) for the quantities Cj(s) in 



terms of the sums fjk, we rewrite the hyperbolic trigono- 
metric functions resulting from the sums in terms of 
exponential functions, and expand the resulting expres- 
sions in powers of exp(— s/2) or exp(— s/4). We also let 
c/L = 1 at this point in the calculation. The time de- 
pendence of the system is recovered by a term-by-term 
inverse Laplace transform of the expansion. The step 
function turn-on of the resulting time dependence arises 
because of the factors exp(— ns/4) in the expansion. The 
lowest order terms in our expansions of the Laplace trans- 
forms are given here: 



2 2e- s / 2 7l [4s 2 + 4s 72 - 4S 2 + tj| - 4i6(2s + l2 j\ 

Cl{S) _ 2s + 7l + {2 S + ll ) 2 [As 2 + 2 S { l2 -i45)-i2 l2 8-48 2 +ui 2 R \ + "' ( ' 

- n 4e-*/ 4 V7i72-(s ~ iS) 

C2l ' SJ " (2s + 71 ) [4s 2 + 2572 - 4(52 +^4-^(45 + 72 )] 



2e _s/2 V7i73(4s 2 - 45 2 +u%- i8s6) 
(2s + 7l )(2s + 73) [4s 2 + 2s 72 - 4<5 2 + w 2 - i25(4s + 72) 



i2e s / 4 ^/7i72o;fl 

d(s) = (2s + 7l ) [4s 2 + 2s 72 - 4<5 2 + co 2 R - i2<$(4s + 72)] + " ' ( 13) 
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r n »2 glfc iAe s/4 ^Wfy(s ~ iS)g 2k 

(s + ikTr)(2s + 71) (s + iibr)(2s + 7l ) [4s 2 + 2s 72 - 4<5 2 +w| - i2(5(4s + 72 )] " 



The inverse Laplace transform of these expressions gives Eqs. (£3| -|27 



APPENDIX B: COEFFICIENTS IN ANALYTICAL SOLUTION 

The coefficients that appear in Eqs. ( p3|j27| ) are given by the following expressions: 

7i72(-7i 2 + 4 l 5 2 +^|-z47 1 J) 
[ 7l 2 -7i72-45 2 +o;| + z2 ( 5(2 7l -72)] 1 ' 

7i [7i 2 ~ 7i72 - 45 2 + tj| + ^(2 7l - 72 )] [ 7l 2 - 2 7l72 - A5 2 + uj 2 r + i4<S(7i - 72 )] 
2 2[7 1 2 -7 1 7 2 -4<5 2 +i4 + i <5(2 7l -72)] 1 j 



^ = 2 7l 7 2 (-7 2 + 77 2 -4t4) (B3) 



V"f2 - (27i ~ 72 + y/ii~ ^R + i^y 



T, = 2 7 i72( 7 2 + V72 2 - 4^) (B4) 



a/72 ~ 4w fl. (- 2 7i + 72 + VtI - - $4$) ' 



77172(71 + »2<5) 
H - 7172 -4^ 2 +w 2 +?;25(2 7 i - 72 ) 



^ - .,2_„,.„,__/|X2 . ,.,2 , .-OJt/O.,. ~ ( B5 ) 



V7T72 (72 - x/t! - 4w ?? ) 

T§ = - ^ y (B6) 

\/7 2 2 - 4c4(2 7 i - 72 + \/7 2 2 ~ 4W 2 , + i45) 

V7172 (72 + V72 - ^r) 

Tn = , ^ ' (B7) 

V7 2 2 - 4^^(-2 7l + 72 + \/7 2 2 " 4^4 - OS) 

~ 4V7i73'(7i 2 ~ 4^ 2 + u R + 

8 — 7 ' \~? r (t>o) 

(73 - 7i) (2 7l - 72 + x/tI - 4w| + i46J ^-271 + 72 + vt! ~ 4t 4 ~ * 4<5 J 

2\/7i7372 (72 - \Jl\ - 4w| J 

^9 = V r-p * ^ (B9) 

\/7 2 2 - 4 ^ (271 - 72 + \/7 2 2 - + *4<5J (- 72 + 2 73 + ^72 - Auj' r + iA5j 

2V2717372 (72 + a/tI ~ 4 <4 ) 

^10 = ; V * ; r (BIO) 

7 7 2 - 4c4 (-2 7l + 72 + 7 7 2 - 4u£ - i45J (73 - 2 73 + 7 7 2 - 4w 2 - i45) 

r 4^/7173 [73 + i(2S + w R )] [73 + i{25 - u R )] 

11 = 7 \— 7 ^=^= r i^ 11 ) 

(73 " 7i) (-72 + 2 73 + v/7 2 2 - 4u R + z4(5j (j 2 - 2 73 + - 4W 2 , - i4£j 



V 7 l 7 2^fl 

°7? - 7i72 - 4<5 2 + uo 2 R + i2(5(2 7 i - 72 ) 



^ = - ^' 7 U ^-. - (B12) 



-r • 2^/ 7l72 w fl 

^13 = « 7 x- (B13) 

Vli - 4[ 4 ( 2 7i - 72 + x/72 2 - M + * 4 <5J 

= 7 2 ^"". c (B14) 



~ 4 <4 (~2 7 i + 72 + VtF 7 M - «45) 



^15 = 

•7"16 = 

Tvt = 

3~18 = 



19 



7i — i2/c7r 



8-^/7172(71 + * 2 <S) 



( 7l - z2fc7r) (2 7l - 7 2 + - 4J| + \A&} (-27! + 72 + ^7| - M - i4<5 

I6 1 /7i7 2 (<5 + far) 



( 7l -i2far) -72 + Vt| - 4w| + i4(<5 + far) 72 + x/tI - 4w| - i4(<5 + far) 
4^7172 (-72 + V72 - 4w h) 



V7l - 4a£ 


-72 + ^72 - 4^4 4 
4^7172 


- i4(8 + far) 
(72 + \/72 ^ 


(27i 


- 72 + VH ~ 4uj r + i46 ) 


^72 - 4^4 


72 + V72 - 4w| - 1 


4(J + far) 


(-271 - 


f 72 + ^72 - 4w£ - z4<5) 
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APPENDIX C: CLASSICAL SCATTERING OF 
EXPONENTIAL PULSES 

A general pulse E(z,t) may be written in terms of its 
Fourier transform as 

1 f°° 

E(z,t) = -= A{u))S k ^'- ur ^ du>, (CI) 

V27T J-00 

where 

1 r°° 

A{uj) = ^ E(0,ty ut dt. (C2) 

V27T J-00 

We consider an incident pulse like that of Eq. ( ^p| ) which 
for convenience we consider to arrive at the origin z = 



at t = 0. The Fourier transform of this pulse is 



A{u) = 1 1 (C3) 

V27T [71/2 - 1{W - Wl)J 



The effect of weak scattering on a Fourier component of 
a plane wave is contained in Eq. (Q), and we construct 
the transmitted pulse to the right of the scattering region 
from the original Fourier components appropriately mod- 
ified according to this equation. Weak scattering means 
that the dimensionless parameter N Az\d\ 2 u> / (heoc^) is 
small, and for convenience we label this parameter /. 
The transmitted field is thus 



E t {z,t) 



2tt J- 



A(oj)e 



iuj{z /c—t) 



8 l2 /2 

+ IJ uj 2 r /4- 52-^72/2 



1 



C_ 

2vr ./_«, [7i/2- 
Ei(z,t) + 



Ju(z/c~t) 



1 + if 



(u - ^2)72/2 



u R /A — - uj 2 ) 2 - i(w - ^2)72/2 



47T 



[71/2 - i(u) - uji)] 



Jui(z/c—t) 



(UJ - ljj 2 ) 



w R /A — (cj — W2) 2 — i(w — ^2)72/2, 



dw. (C4) 



The integral in the last line above gives the scattered field 
E s (z,t). We evaluate the integral using contour integra- 
tion, and note that all three poles are in the lower half 
of the complex plane. For (z/c—t) > the integration 



along the real axis is closed in the upper-half plane, giv- 
ing a result of zero; for (z/c — t) < the contour is closed 
in the lower half plane encircling the poles. Evaluating 
this integral gives 
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E s (z,t) = Cf l2 e 



(71 +i2<5)e-^(*- z / c ) 



[ 7 2 - 7! 72 - AS 2 



- 126(2^1 - 72 



(72 - v ^ 2 ~4^|) e K^-V^- 4 ^)/ 4 +^( t - 2 / c ) 
)] V72 ~ 4w fi ( 2 7i " 72 + ^72 ~ 4t 4 + i4(5 ) 



+ 



( 7 2 + ^34^) e [-(72+V^4^)/4+^(t-»/c) 

a/72 - 4t 4 ( 2 7i - 72 - ^72 - 4t 4 + i4(5 ) 



(C5) 



This is identical in form to the terms describing the 
scattered field in the quantum mechanical expression 
Eq. ©. (The factor e~ lUlt does not appear in Eq. (§§) 
because of the zero chosen for the energy scale in the 
quantum calculations.) 
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